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Abstract 

The measure of the bullwhip effect, a phenomenon in which demand vari- 
ability increases as one moves up the supply chain, is a major issue in Supply 
Chain Management. Although it is simply defined (it is the ratio of the 
unconditional variance of the order process to that of the demand process), 
explicit formulas are difficult to obtain. In this paper we investigate the theo- 
retical and practical issues of Zhang [Manufacturing and Services Operations 
Management 6-2 (2004b) 195] with the purpose of quantifying the bullwhip 
effect. Considering a two-stage supply chain, the bullwhip effect is measured 
for an ARMA(p,q) demand process admitting an infinite moving average rep- 
resentation. As particular cases of this time series model, the AR(p), MA(q), 
ARMA(1,1), AR(1) and AR(2) are discussed. For some of them, explicit for- 
mulas are obtained. We show that for certain types of demand processes, the 
use of the optimal forecasting procedure that minimizes the mean squared 
forecasting error leads to significant reduction in the safety stock level. This 
highlights the potential economic benefits resulting from the use of this time 
series analysis. Finally, an R function called SCperf is programmed to cal- 
culate the bullwhip effect and other supply chain performance variables. It 
leads to a simple but powerful tool which could benefit both managers and 
researchers. 
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1. Introduction 



In recent years, companies in various industries have been able to signifi- 
cantly improve their inventory management processes through the integration 
of information technology into their forecasting and replenishment systems, 
and by shar ing demand-related information with their supply chain partners, 



Avivl ( ]2003[ ). However, despite the benefits resulting from the implementation 



of the above practices, inefficiencies still persist and are reflected in related 
costs. 

The bullwhip effect, defined as the increase in variability along the sup- 
ply chain, is a frequent and expensive phenomenon identified as a key driver 
of inefficiencies associated with Supply Chain Management (SCM). It dis- 
torts the demand signals, which causes instability in the supply chain, and 
inc reases the cost o f supplying end-customer demand. 



Forrester! (119581 ) was the first to popularize this phenomenon. Inspired 



by Forrester's work, several researchers have studied the bullwhip effect. 



Stermanl ( 119891 ) used the Beer Game, the most popular simulation of a simple 



production and distribution system, to demonstrate that the bullwhip effect 
is a significant problem with important managerial consequences. It results 
in unnecessary costs in supply chains such as inefficient use of production, 
distribution and storage capacity, recruitmen t and tra i ning; costs, increased 



i nvento ry and poor customer service levels ( iMettersI (119971 ) and iLee et al 
(Il997bh). , , 



Lee et al.l (jl997al Jbl) identified four main causes of the bullwhip effect: 



demand forecasting, order batching, price fluctuation and supply shortages. 
Of these, demand forecasting is recognised as one of the most important since 
the inventory system is directly affected by the forecasting technique chosen. 
Three popular forecasting methods are commonly used: the Minimum Mean 
Squared Error (MMSE), Moving Average (MA) and Exponential Smoothing 
(ESI ., , 



Chen et al.l (j2000al ) quantify the bullwhip effect considering the MA fore- 
cast method for a simple two-stage supply chain and a first-order autoregres- 
sive demand process, AR(1). The authors show that the bullwhip effect is in 
part due to the effects of demand forecasting. Therefore, given complete ac- 
cess to customer demand information for each stage of the supply chain, the 
bullwhip effect can be significantly reduced. However, they also show that 
the bullwhip effect will exist even when demand information is shared by all 
stages of the supply chain and all stages use the same forecasting technique 



2 



and the same inventory policy. In similar work IChen et al.l (l2000bl ) quantify 
the bullwhip effect considering this time the ES forecast and two different de- 
mand processes: AR(1) demand process and a demand process with a linear 
trend. In both works, the authors recognize an important limitation of their 
results: the models considers only non-optimal forecasting methods. The 
authors justify this limitation saying that ES and MA are commonly used 
in practice. Users are in general less familiar and less satisfied with more 
so phistica t ed met hods like time series techniques. 



Zhangl (j2004al ) investigates the impact of MMSE, MA and ES forecasting 



methods on the bullwhip effect for a simple inventory system in which AR(1) 
demand process describes the customer demand and an Order-Up- To (OUT) 
inventory policy is used. The study shows that different forecasting methods 
lead to bullwhip effect measures with distinct properties in relation to lead- 
time and the underlying parameters of the demand process. The author 
shows that MMSE forecasting method leads to the lowest inventory cost. 
This result is not surprising since MMSE method is optimal when the demand 
model is known to be an AR(1) process. On the other hand, if the demand 
structure is not well known, the MA or ES method may perform better than 
the MMSE method because they are more flexible. 

Another aspect studied in relation of the bullwhip effect is the demand 
process. A variety of time-series demand models have appeared in the liter- 
ature of inventory control and SCM. By far, the AR(1) process is the most 



frequen tly adopted demand model to study the b ullwhip effect (IChen et al. 



feOQOaU bh. iLee et al.l fll997al Jbh and|ZhanJ <j2004ah ). Rece nt works use more 
sophi sticated time series models like ARM A an d ARIMA (|Box and Jenkins! . 
1970l ) to have more realistic demand models. Luong and Phien (2007) use 



an AR(2) and a general AR(p) model; iDuc et al.l (120081 ) use an ARMA(1,1) 
model. In all these models an analytical derivation of the bullwhip effect 
measure is presented and the effects of the autoregressive coefficient on the 
bul lwhip effect is investigated. 



Zhangl ( 12004bl ) uses an ARMA(p,q) model to study the demand evolu- 
tion in supply chains. The author shows that the order history preserves 
the autoregressive structure of the demand. Zhang's work identifies an im- 
portant application of this result relating to the quantification of the bull- 
whip effect. In this paper, inspired by Zhang's work, we study the the- 
oretical and practical issues in order to measure the bullwhip effect for a 
generalized demand process. In addition, we programmed a function in R 
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( R Development Core Team . 201(j| ). called SCperf 1 , which implements the 



bullwhip effect and others supply chain performance variables. It is well 
known that measuring the bullwhip effect is difficult in practice but the 
SCperf function overcomes this problem thanks to the help of an R func- 
tion (ARMAtoMA) which converts an ARMA process into an infinite moving 
average process. As far as practical applications are concerned, the economic 
implications of this phenomenon on the inventory cost have been considered. 

Our contributions to this subject can be described as follows: first, this 
study hopes to improve the understanding of time series techniques. Second, 
we show that for certain types of demand processes the use of the opti- 
mal forecasting procedure that minimizes the mean squared forecasting error 
leads to significant reduction in the safety stock level. This highlights the 
potential economic benefits resulting from the use of this time series analysis. 
Finally, the SCperf function leads to a simple but powerful tool which can 
be helpful for the study of this phenomenom and other supply chain research 
problems. 

The structure of our paper is as follows. The next section presents 
the inventory model. Section [3] presents a general ARMA(p,q) case with 
ARMA(1,1), MA(q), AR(p), AR(1) and AR(2) as particular cases. Next the 
economic implications are shown. The final section summarizes the main 
results of the research. 

2. Inventory model 

In this paper we consider a simple supply chain model for a single item 
and an OUT inventory policy in which the retailer determines a target level 
or OUT level and, for every review period, places a n order suffic i ent to b ring 



the inventory position back to this level. As did I Chen et al.l (|2000bJ), we 
consider that the ordered quantity made in period t is received at the start 
of period t + L where L is defined to be a fixed lead time plus the review 
period, i.e., L is the lead time plus 1. For instance, in the case of zero lead 
time, L = 1. Shortages are back-ordered and no fixed ordering cost exists. 
In the remainder of the paper L will call the lead time. This choise is made 
for sake of brevity, and should not create confusion. 

The sequence of events during a replenishment cycle for each period t can 
be described as follows: the retailer receives orders made L periods ago; the 



See the supplementary material 
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demand d t is observed and satisfied; the retailer observes the new inventory 
level and finally places an order O t to the supplier. As a consequence of this 
sequence of events, the ordered quantity can be written as: 

O t = S t - St-i + dt, (1) 

where St represents the OUT level in period t, i.e., the inventory position at 
the beginning of period t. Note that in the above expression, we have implic- 
itly assumed that the order quantity can be negative, i.e., returning items 
are allowed at no costs. This unpleasant feature is needed for tractability. 
However, the free- return assumption becomes negligible when the demand 
m ean is suffi c iently larg e. Further deta il abou t this assumption can be found 



in 



Lee et al.l ((20001) and IChen and Led (120091 ) . 



Under the OUT policy, the OUT level St can be estimated from the 
observed demand as: 

St = Dt + zd}, (2) 

where D\ = Xlr=i dt+r is an estimate of the mean demand over L periods 
after period t, z is the safety factor which is a fixed constant chosen to meet 

a required service level and of = JvariDf - £%) is an estimate of the 
standard deviation of L periods forecast error. An OUT policy of this form 
is optimal when the demand came from a normal distribution and there is 
no set up or fixed o r der cos t. 



As I Chen et al.l (l2000bl ) mention, if the retailer follows an OUT policy 
of the form S t = D L + za L , where D L is the known mean and a L is the 
standard deviation of the demand over L periods, then the OUT level in any 
period is constant and, consequently, the order is equal to the last observed 
demand. Therefore, there is no bullwhip effect. However, these values are, in 
general, unknown and the retailer must estimate them using some forecasting 
technique. Note that the introduction of forecasting values in the calculation 
of St is one of the main causes for the variability increase along the supply 
chain or, in other words, the bullwhip effect. 

The demand forecast is performed here by using the MMSE method. 
It was shown that, for an ARMA process, the MMSE forecast for period 
t + t is the conditional mean given the observed information 2 . Let F t = 



2 Box and Jenkins, 1970, pp.128. 
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{d t ,d t -i, ....} be the information set which represents all the information 
available until period t. Hence, the demand forecast for r periods ahead is 
given by E(d t+T \F t ). 

In order to quantify the bullwhip effect we combine and (J2J) to rewrite 
the order quantity as: 

O t = (Df - + A° L t - SjLi) + d t . (3) 

We show later in the paper (see Lemma IA.1I) that the standard deviation 
of lead-time forecast error remains constant over time for an ARMA(p,q) 
demand process. Hence, a\ = a\_ x and the order quantity given in (j3J) 
becomes 

o t = {b h t - b L t _ x ) + d t . (4) 

Let M be the measure for the bullwhip effect. Since M can be obtained 
from the ratio between the unconditional variance of the order process to 
that of the demand process, we have 

m = ( 5) 

Var(d t ) y ' 

Note that M is calculated by using the variances from both side of Equation 
PJ. The fact that M = 1 means that there is no variance amplification, while 
M > 1 means that the bullwhip effect is present. On the other hand, M < 1 
means that the orders are smoothed if compared with the demand. The last 
case is less common since it is unlikely to have a situation where stages up 
the supply chain have a better representation of the customer demand than 
the first stage (i.e., the retailer). 

In what follows, the corresponding bullwhip effect measure is derived for a 
general ARMA(p,q) demand process and some particular cases are discussed. 
Since the calculation is complex, we cannot always express this measure in a 
closed form. In this context, the SCperf function was developed to overcome 
this computational difficulty. 



6 



3. ARMA(p,q) case 



The demand process, d t , seen by the retailer, is described by a stationary 
ARMA(p,q) process as follows 3 : 

d t = /i + H h 4>pd t -p + e t + 9 x e t -x H h 6 g e t - q , (6) 

where /i is a nonnegative constant, e 4 is i.i.d. normally distributed, with 
mean zero and variance of, p is the autoregressive order of the process, g is 
the moving average order of the process, <pj is the autoregressive coefficient, 
and 8j denotes the moving average coefficient. It is often useful to express 
fl6]) in terms of the lag operator, B, where B k d t = dt-k- In order to do so, let 

<p(B) = 1 - faB <p p B p and 6(B) = 1 + 9 X B H h 6 q B q . Hence, the 

demand process in (jfJJ) can be expressed as: 

<f>(B)d t = {i + 9(B)e t , 

where 4>(B) and 9(B) are known as the autoregressive and the moving average 
polynomials in the lag operator of degree p and q. If we substitute the lag 
operator by a constant z, we get the characteristic equations: 

<p(z) = 1 - - (f) 2 z 2 4> p z p 

and 

9(z) = l + 9 lZ + 9 2 z 2 + ■■■ + 9 q z q . 

The process is called the autoregressive process of order p, AR(p), if 9(z) = 1 
and a moving average process of order q, MA(q), if <p(z) = 1. We assume 
that the process described in ([6]) is invertible and covariance stationary, i.e., 
the roots of the equations 9(z) = and (j>(z) = must be outside the unit 
circle. To avoid the problem of parameter redundancy, it is assumed that the 
two characteristic equations share no common roots. 

It is important to note that the constant z in the above equations is 
different from the constant used to define the safety factor. We have chosen 



3 Our representation differs from some works where the MA model is written with 
negative coefficients, i.e., dt — fJ, + 4>id t -i + ■ • • + 4> p dt- p + e t — #i£t-i — • ■ ■ — 9 q e t - q - 
We chose this representation to be in accordance with the R software which was used to 
implement the bullwhip effect. 
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this notation to be in accordance with time series notation and we hope 
that this will not cause any future confusion. Using stationarity and taken 
expectations in fl6]) directly it can be found that the mean of ARMA(p,q) 
demand process is defined by 

1-01 <Pp 

It is known from time series theory that a stationary ARMA(p,q) demand 
process under the above conditions can be written as an infinite moving 
average process of its errors, MA(oo), that is, 

<k = fi d + ££= ^e H) (8) 
where fid is defined as in Equation (0) and the sequence {ipj} in (jHJ) is de- 
termined by the relation tp(z) = Y^Lo^j^ = ^?f)> or equivalently by the 
identity 

(ip + ^ 1 z+^ 2 z 2 + - ■■)(l-<p l z-<p 2 z 2 <p p z p ) = (l+6 1 z+d 2 z 1 + - ■ - + 9 q z q ). 

Equating coefficients of z\ j = 0, 1, we find that 

p 

= X] Mj-k + &j for 3 > 1> ( 9 ) 

k=l 

where 9 = 1, 9j = for j > q, and ipj = for j < 0. Note that equation (jHJ) 
is a recursive equation. Therefore, the -^-weights satisfy the homogeneous 
difference equation given by 

p 

$i ~ Mj-k = °> 3 > max(p, q + 1), (10) 

k=l 

with initial conditions given by equation (Q. From homogeneous difference 
equation theory the general solution for equation ( TTOT) can be read off directly 
as: 

Ipj = cxz[ ] H hc r z~ J , (11) 

where distinct roots of the polynomial <p(z) and Ck, for fc = 

l,2,...,p are constants which depend on the initial conditions. 4 Now, from 



4 In the case of the repeated root, the solution is different. See IShumwav and Stoffer 
( 20061 ) for a brief and heuristic account of the t opic. Fo r deta ils about homogeneous 
difference equation theory the reader is referred to iMickensI (|l987l ) . 
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equation (JSJ), the variance of the demand process can be expressed as: 



3=0 



(12) 



It is important to note that the MA(oo) representation depends on an 
infinite number of parameters and, conse quently , it is n ot directly useful in 
practical applications. On the other hand, IZhang (2004b)), using the MA(oo) 
representation, shows a property, called by the author ARMA-in-ARMA-out 
(AIAO), which reveals that the order history preserves the autoregressive 
structure of the demand and transforms its moving average structure accord- 
ing to a simple algorithm 5 . As the author remarks, the practical value of 
the AIAO property lies in its ability to make simpler the measuring of the 
bullwhip effect. 

Proposition 1. (Zhand , 200$) The retailer's demand process can be rep- 
resented by an MA(oo) process with respect to the retailer's full information 
shocks e t , as in equation |S)). Hence, the retailer's order O t to its supplier is 
given by: 



3=0 3=1 

where the x/jj = for j < 0, ipo = 1, and ipj = Y^k=i ^k^j-k + &j for j > 1 



Proof. See lZhand fl2004bl ). 



(13) 



□ 



Proposition 2. For a stationary ARMA(p,q) demand process, the measure 
for the bullwhip effect is defined by: 



M = 1 + 



j=i+l 



(14) 



where the ipj = for j < 0, ip = 1, and if)j = Y^k=i ^k^Pj-k + 6j for j > 1. 

Proof. Taking the variance of the order quantity, Equation ( JIB"]) , we have 

Var(O t ) = o-^ZU ^) 2 +^ 2 £r=i = ^T=o ^I+ 2 Eto E jUn Mil 
We complete the proof by substituting this result and (fT2j) in ([5]). □ 



Zhang 2004b, pp. 197 
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Proposition 3. The bullwhip effect increases when the lead-time L increases 
if and only if ip L+1 Y,j=o > °- 



Proof. From equation (JHJ), it is straightforward to see that the bullwhip 
effect exists, i.e., M > 1, if and only if Ei^=o Ej=j+i ^i^j > 0- Let g{L) = 
Ef=o Ej=i+i and A #( L ) = 9(L + 1) - g(L). Then 

A #(£) = E^Ej^i^i-ELjEiU+i^i = ^o(E J L =iVj-Ei=iV'i)+ 

h ^L-i(Et=i ^3 ~ ^l) + = ipL+i EjU V'i- Hence, A#(L) > if 

and only if ^l+i E^=o > ^- Hence, is a non-decreasing function of 
the lead-time L if and only if ipL+i Ej=o > 0- '— ' 

5.1. ARM A( 1,1) case 

The stationary ARMA(1,1) demand process is described as follow: 

df = /i + 4>d t -i + e t + 9e t -i. (15) 

Stationarity and invertible conditions impose \(f)\ < 1 and \8\ < 1. It can be 
shown that the mean and variance of the demand process are fid = and 

a 2 _ ( 1+e i 2 +^f ) g ? } respectively. 

Proposition 4. For a stationary ARMA(1,1) demand process the measure 
for the bullwhip effect is defined by: 

M < L - *• »> = 1 + { i%Z>tm I 1 - ^ t+I + 9 ^ - *">] • < 16 > 

Proof. Since the AR polynomial associated with (ITST) is </>(z) = 1 — <pz, and 
its root, say z±, is z\ = then the general solution for the ^-weights can 
be written directly from equation (ITT)) as ipj = c(j> 1 . From (|9]) we find that 
the initial conditions are ipo — 1 and ipi = <ft + 9, which combining with the 
general solution, results in c = (<fr + 6)/<f). Hence, ipj = (<f) + O)^^ 1 for j > 1. 
Since we know iftj, we can rewrite the follow relations as: 

L L ILL 
i=0 j=i+l j=l i=l j=i+l 

= + wT + (i-*)(i-^) 

(0 + #)(i-<^ 



0)(1-0 2 ) L ^ 
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-0.5 



0.0 



Figure 1: Bullwhip generated with ARMA(1,1) demand process when L=l 



and 



E 

3=0 



i + e 2 + 2<pe 



Substituting the two above results in equation (flip we complete the proof. □ 
Using a generalized formula for the va riance ratio, we get a similar expression 
to that obtained by iDuc et al.l ( 120081 ). There are two other results found by 
the above authors which are easily verified. 



Proposition 5. The bullwhip effect exists, i.e, M(L,<f>,0) > 1, if and only 
if (j) + 9 > 0. 



Proof. IDuc et al.ll2008l . pp. 248-249. 



□ 
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Proposition 6. The bullwhip effect, measured by M(L,<j),9), has the follow- 
ing properties. 

(a) If <f) > 0, the bullwhip effect increases as L increases. 

(b) If —9 < (f) < and L is an odd number, the larger L is, the smaller 
the bullwhip effect is. 

(c) If —9 < (p < and L is an even number, the larger L is, the larger 
the bullwhip effect is. 



Proof. IDuc et al.ll2008l . pp. 249. □ 



In conclusion the bullwhip effect occurs only when the sum of the AR 
parameter and the MA parameter is larger than zero ( See Figure [1]) and it 
does not always increase when the lead time L increases. In fact, if (p + 9 > 
and > the bullwhip effect increases when the lead-time increase. However, 
if —9 < (f) < and L is an odd number, the bullwhip effect becomes smaller 
as L becomes larger; if —9 < (j) < and L is an even number, the bullwhip 
effect becomes larger as L becomes larger. Figure [2] represents situations 
where these facts are observed. 

3.2. MA(q) case 

The MA(q) demand process can be written as 



g 

d t = n + J2 Ofr-j =fJi + (l + 6 1 B + --- + 9 q B q )e t = fi + 9(B)e t . 

Since 9(B) is finite, no restrictions on the MA parameters are needed to 
ensure stationarity. Considering q — > oo the infinite MA representation is 
written as: 

oo 

d t = /id + y^j-et-j, 

3=0 

where ipj = 9j for j = 0, 1, .., q and ipj = for j > q. It can be easily seen 
that fid = A 4 an d v\ = (1 + 0? + " " ' + $g)°f • Since the above demand process is 
i.i.d. the OUT level, St, is constant across all periods. Hence, from Equation 
([!]), O t = d t , consequently, the bullwhip ratio equals one. 
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theta=-0.6 



theta=0.1 





3.3. AR(p) case 

The stationary AR(p) demand process is described as follow: 

d t = fi + 4>id t _i H h 4>pd t ^ p + e t 

Assume that the AR parameters are such that {d t } is stationary. It is 
straightforward to verify that the MA(oo) representation is 

d t = ftd + i>(B)e t , 

where fid is defined as in (JTj) and ip{B) = The -^-weights in the 

MA(oc) representation of d t are found directly from (fTTj) and it can be shown 
that the constants are expressed by: 

z p-i 
13 



wher e the constants terms q sum to the unity, cH — - + c p = 1, see lHamilton 



1994 pp. 33-36, for details. 



3.4. AR(1) case 

The stationary AR(1) demand process is described as follows: 

d t = 11 + 0d t _i + e t . (18) 

Stationarity condition imposes \4>\ < 1. Using stationarity it can be shown 
that the mean and the variance of the process are fid - 
respectively. 



and aj = ^ 



Proposition 7. For a stationary AR(1 ) demand process the measure for the 
bullwhip effect is defined by: 



M(L, 



20(1 - (j) L )(l - (f) L+1 ) 



(19) 



Proof. As in the ARMA(1,1) case, the AR polynomial associated with (ITS]) 
is (f)(z) = 1 — 4>z, and the root, say, z±, is z\ = <p~ l . Using (fTTl) the general 
solution is ipj = c{z\)~i = c<f>\ with ip = 1 and = as initial conditions. 
Combining the general solution with the initial conditions we find ipj = <fP . 
Since ipj = 



Equation f jMj) can be expressed as: 



M(L, 



2EtoE 



j=i+l 



where 



EE w 

i=0 j=i+l 



L L-i-l 

EE> 

i=0 k=0 



Liik+i+1 



E 

i=0 



1-0 2 1 - 

(1 -0 L )(1 -0 L+1 )" 



(20) 



and Y%Lo 1 
the proof. 



l-4> 2 



. Substituting the two above results in (120]) complete 

□ 
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T 

-1.0 -0.5 0.0 0.5 1.0 

phi 

Figure 3: Relationship between the bullwhip effect and demand autocorrelation 

Proposition 8. For a stationary AR(1) demand process the bullwhip effect, 
measured by Equation [Tty), has the following properties: 

(a) The bullwhip effect exists, i.e, M(L, 0) > 1, if and only if(f)>0. 

(b) For cf) > 0, a longer lead-time leads to a more significant bullwhip 
effect. 

PROOF. Since 1 - > 0, 1 - L > and 1 - L+1 > for |0| < 1, 
it is straightforward to see that M(L, 0) > 1, if and only if > 0. Let 
/(L,0) = 0(1 - L )(1 - i+1 ) and Af(L) = f(L + 1,0) - /(L,0). Then, 
A/(L) = (1 - 2 )(1 - i+1 )0 L+1 . It can be easily seen that A/(L) is an 
increasing function with respect to L since > 0. Hence, the bullwhip effect, 
i.e, M(L,0), increases as L increases since > 0. □ 

Figure [3] depicts how the bullwhip effect generated by AR(1) demand process 
increases for different lead-time values, L = 1, 6. We can observe that the 
increase of the lead-time has a strong impact on the bullwhip effect when 
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> 0.5 and a less significant one wh en is positive and near zero and one. 



Therefore, as it was already noted by [Zhang (l2004al ). reduction on the lead 



time can reduce the bullwhip effect if the demand autocorrelation is positive 
and away from zero and unity in the case of AR(1) demand process. 

3.5. AR(2) case 

The stationary AR(2) demand process satisfies: 

dt = \i + 4>i<k-i + 4>2<k-2 + e* (21) 

In the AR(2) case, stationarity implies that the roots of <j)(z) = lie outside 
the unit circle or, equivalently, the parameters 0i and 02 must lie in the 
triangular region restricted by 0i + 02 < 1, 02 — 0i < 1 and |02| < 1- It can 
be shown that for a stationary AR(2) demand process the mean and variance 
of the demand are ^ an d (1+ ^y_ 2 ^-^p respectively. 

Proposition 9. Let z\ and zi be the solutions for the characteristic equation 
defined by the AR(2) process. For a stationary AR(2) demand process the 
■ip-weights are defined by: 

l+j _ J-+3 

if). = ^ z -l 

ZiZ 2 (z 2 - Z-l) 

Proof. From Equation ( ITOj) . the general solution for ^-weights for an 
AR(2) process is described by: 

^ = Cl ( Zl )- j + c 2 (z 2 )-i (22) 

where 

= - h ± 

202 

and 

202 

are the solutions for the characteristic equation 1 — <piz — <p 2 z 2 = 0. On the 
other hand, from Equation (IT7|) . the values of the constants are given by: 



ci = _/ zi (25) 
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and ^ 

C2 = - J 2 _ -x (26) 

z x z 2 

Finally by replacing (|25]h pi]), p5|) and (p6) in (122]) we find the result. □ 

Note that the solution for the ^-weights are a function of the roots of the 
AR polynomial. In the AR(2) case, the roots can be real if <p\ + 40 2 > 0, 
or complex if <fi\ + 40 2 < 0. In both cases, from a computational point of 
view, the solution for the V^-weights can be found and, therefore, we can get 
a measure for the bullwhip effect. Since an explicit form for the measure for 
the bullwhip effect is difficult to obtain, we investigated the relation of the 
autoregressive coefficients and lead-time by num erical experimentation. For 



an analytical derivation the reader is referred to iLuong and Phienl (120071 ) . 

When 0i < 0, the bullwhip effect does not exist for 2 < and for 2 > 0, 
02 — 0i < 1- On the other hand, when <pi > the bullwhip effect always 
exists for 2 > 0, 0i + 2 < 1 and for 2 < 0, 0i + 2 < 1. The pattern 
shown when the lead-time is equal to one does not seem to be the same when 
the lead-time increases. Using the function SCperf , it can be verified that 
the there is no bullwhip effect when X < and 2 < and always does 
when 0i > 0, 2 > and 0i + 2 < 1. In the last case, we observe that the 
bullwhip effect increases when the lead-time L increases, see Table [TJ 

Table [T] also shows that there is no clear relation between the autoregres- 
sive parameters and the bullwhip effect when they have different signs. In 
these situations the bullwhip effect may or may not exist depending on the 
values of 0i, 2 and L, and it does not always increase when lead-time in 



crease s. These remarks confirm the results pointed out by ILuong and Phien 



( 120071 ) . 

In conclusion, when both first-order and second-order AR parameters are 
positive, the bullwhip effect exists and it increases as lead-time goes up. 
However, when the AR parameters have different signs the behaviour of the 
bullwhip effect is not clear. The bullwhip effect does not always exist and 
it is not always correct that the bullwhip effect necessarily increases when 
lead-time increases. 



4. Economic implications 

An important economic application of the use of time series methods can 
be seen in the safety stock level, which is the amount of inventory that the 
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Table 1: Bullwhip effect generated for different AR(2) demand pro- 
cess.* 



L AR(c(-0.2,0.7)) AR(c(0.6,-0.4)) AR(c(0.7,0.2)) 



1 


0.886667 


1.822857 


1.315000 


2 


1.222133 


1.735086 


1.842850 


3 


0.970805 


1.170277 


2.512887 


4 


1.379174 


0.917179 


3.291280 


5 


1.051166 


0.949074 


4.141105 


(3 


1.450366 


1.060235 


5.035836 


7 


1.097494 


1.117111 


5.953552 


8 


1.464249 


1.103809 


6.877221 


9 


1.117408 


1.072652 


7.793541 


10 


1.447477 


1.059437 


8.692330 



* SL=0.95 



retailer needs to keep in order to protect himself against deviations from 
average demand during lead time. 

Let SS = zadVL and SSLT = zo\ be two safety stock measures. The 
former is traditionally used in some operational research manuals and it is 
based on the standard deviation of the demand over L periods, the latter is 
the safety stock as defined in (j2J) and it is based on the standard deviation 
of L periods forecast error. 



Chen et al.ll2000bl pp. 271, pointed out that SSLT will be greater than 



SS, i.e., using time series analysis, the retailer will hold more safety stock to 
achieve the same service level. According to the authors this is because SS 
captures only the uncertainty due to the random error e and SSLT captures 
this uncertainty plus the uncertainty due to the fact that the mean demand 
D\ is estimated by Df, in our case using the MMSE forecasting method. We 
show by numerical experiments that for some special cases SSLT is lower 
than SS regarding lead-time and service level. 

Using the SCperf function, it was verified that for ARMA and AR cases, 
high values on AR parameters and small values of lead-time result in lower 
SSLT . However, in general, there is a lead-time value for which this situation 
is reversed. Table [2] shows the safety stock levels SS and SSLT generated 
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Table 2: Bullwhip, SS and SSLT generated 
by ARMA(0.95,0.4) demand process.* 



T 

L 


Bullwhip 


SS 


nriT rp 

SSLT 


1 


1.13711 


7.299 


1.645 


2 


1.44321 


10.323 


4.201 


3 


1.89270 


12.643 


7.304 


4 


2.46294 


14.598 


10.817 


5 


3.13393 


16.322 


14.652 


6 


3.88802 


17.879 


18.745 


7 


4.70970 


19.312 


23.048 


8 


5.58531 


20.645 


27.522 


9 


6.50289 


21.898 


32.137 


10 


7.45199 


23.082 


36.867 



* SCperf(0.95,0.4,L,0.95) 



by ARMA(0.95, 0.4) demand process and service level equal to 0.95 for ten 
different values of lead-time, L = 1, .., 10. For instance, for L = 2 we have 
SS = 10.3 and SSLT = 4.2, a difference of 6 units which represents a saving 
of 59.2% over SS. Note that this difference decreases when the lead-time 
increases until L = 6 where we have SSLT larger than SS. 

It is difficult to know for which value of lead-time SSLT becomes larger 
than SS. In general, it depends on the AR parameters of the demand. For 
negative values of the AR parameters, it occurs for lower values of lead-time. 
Nevertheless, for the AR(2) case the AR parameters present a more complex 
relation with the performance of the SSLT. When the first-order and second- 
order AR parameters are positive, the pattern is the same as the AR and 
ARMA case, that is, SSLT becomes larger than SS for high values of lead- 
time. Moreover, when the first-order and second-order AR parameters have 
different signs, it is difficult to determine when the SSLT is better than SS 
as a measure for the safety stock level. 

Table [2] shows that there is a benefit resulting from the use of SSLT 
instead of SS as a measure for the safety stock level when regarding the lead- 
time. This benefit was verified for special demand processes where the AR 
parameters are high. Moreover, if for those lead-time values where SSLT 
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Table 3: SS and SSLT generated by different demand processes 



Models 


Service Level 


L= 


1 


L= 


-2 


L= 


=3 




SL 


SS 


SSLT 


SS 


SSLT 


SS 


SSLT 




0.90 


5.687 


1.282 


8.043 


3.273 


9.850 


5.691 




0.91 


5.950 


1.341 


8.414 


3.424 


10.305 


5.954 




0.92 


6.235 


1.405 


8.818 


3.588 


10.800 


6.239 




0.93 


6.549 


1.476 


9.262 


3.769 


11.343 


6.553 




0.94 


6.899 


1.555 


9.757 


3.971 


11.950 


6.904 


ARM A(0. 95, OA) 


0.95 


7.299 


1.645 


10.323 


4.201 


12.643 


7.304 




0.96 


7.769 


1.751 


10.987 


4.471 


13.456 


7.774 




0.97 


8.346 


1.881 


11.803 


4.803 


14.456 


8.352 




0.98 


9.114 


2.054 


12.889 


5.245 


15.785 


9.120 




0.99 


10.323 


2.326 


14.599 


5.941 


17.881 


10.330 



is smaller than SS, we consider the service level, it is verified that SSLT is 
always smaller than SS when the service level increases. 

Table |3] presents SSLT and SS generated by the same demand process for 
L = 1,2,3 and ten different values of service level, SL = 0.9, 0.91, 0.99. 
Note that when considering the service level, the difference between SS and 
SSLT increases for larger values of service level differently when lead-time is 
regarded. For instance, for L = 1 and SL = 0.97 we have SS = 8.35 and 
SSLT = 1.88. There is a difference of 6.47 units which represents a saving 
of 77.46% over SS. 

All of these facts suggest that there is a potential benefit resulting from 
the use of time series analysis when regarding the lead-time for some demand 
processes and, in this context, the benefit is even greater when the service 
level is considered. On the other hand, the relationship between the bullwhip 
effect measure and the safety stock level is more complex. Although Table 
[2] shows a positive relation between the bullwhip effect and the safety stock 
level, this relationship is not completely clear as can be seen using the SCperf 
function for the AR(2) case when <fii = —0.2 and 02 = 0.7. 

In conclusion, when inventory cost and service level are of primary concern 
the MMSE forecast should be used since it leads in some cases to lowest safety 
stock level. Although the MMSE forecasting requires more computational 
effort, the SCperf function implements this method in an easy way. 
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5. Summary 

In this paper we quantify the bullwhip effect using Zhang's result for 
a stationary ARMA(p,q) demand process which admits an MA(oo) repre- 
sentation. It is well known that measuring the bullwhip effect is difficult 
in practice. We show that using a generalized form of this measure, the 
computation of this ratio is simplified if compared with traditional recursive 
procedures. In some particular cases we obtain explicit formulas for this 
ratio. 

The SCperf function was programmed in R which implements the bull- 
whip effect. We have evidenced that the use of this function makes possible 
accurate estimations of the bullwhip effect and other supply chain perfor- 
mance variables. We point out that no approximation is required. Moreover, 
we show that for certain types of demand processes the use of MMSE con- 
sidered in the model leads to a significant reduction in the safety stock level 
regarding lead-time and service level. All of these observations highlight the 
potential economic benefits resulting from the use of time series analysis but 
it depends on the underlying demand process. For instance, if we consider 
an ARM A (1,1) demand processes with a high AR parameter, the use of time 
series techniques leads to a significant reduction in the safety stock level but 
this is not the case when a low AR parameter is considered. 

The SCperf function leads to a simple but powerful tool which gives exact 
analytical solutions to a set of supply chain equations, opening up a whole 
new range of research opportunities. Moreover, since the function presented 
in this paper is easy to use, it might be used to complement other managerial 
decision support tools. Finally, the code is given, which makes, together with 
the fact that R is freeware, the whole research reproducible by everyone. It 
may also be modified for specific tasks. 
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Appendix A. 



Lemma 1. For a stationary ARMA(p,q) demand process, the variance of 
forecasting error for the lead-time demand remains constant over time and 
is given by: 



a L t ) 1 = Var(D L t - D 



i + (E^) 2 + --- + (E^-) 2 



(A.l) 



where ipj satisfy (Q and (jTOj) and is given by (fTTl) . 



Proof. Since L> t L = J2t=i dt+ T , A L = Y1t=i dt+r with r — 1, .., L and d t+T = 
E(d t+T \Ft) — + Yl'jLr'ipjtt+T-j, the variance for the lead-time demand 
forecast error is 



L\2 



Var[D^-D^] = Var 

L 



T = l 



Var 



Var 



T=l \j=0 
' L 7—1 

r=l j=0 



By expanding the above double sum and combining the same error terms, it 
follows that: 



L r-1 1 
r=l j=0 j=0 



L-l 
j=0 



The independence of future error terms leads to the variance formula for 
lead-time demand forecast. □ 
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